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Abstract 

This paper presents a comparative analysis of complexes made of one Rare Gas (Rg) atom and 
molecular hydrogen, for all five stable Rg atoms. In particular, the vibrational band origins have 
been calculated, as well as the structural properties of the associated wavefunctions. The study is 
extended to cold Rg-H2 scattering. The molecular systems are studied variationally using a very 
simple, yet effective, trial wavefunction. A large number of Potential Energy Surfaces available from 
the literature is considered. A comparative analysis shows that differences of up to two orders of 
magnitude exist for the zero energy elastic cross sections of the five complexes. Corrections to the 
model have also been considered, showing no significant effect. 
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I. INTRODUCTION 



Over the last decades, complexes made of one rare gas atom and molecular hydrogen 
Rg-H2) have provided an optimum testing ground for ab initio and experimental techniques 
ij. A wealth of data for the Rg-H2 complexes has been collected in a large number of 
very accurate experimental investigations. This includes high-precision measurements 
of the rotational- vibrational transitions (mostly in the microwave and infrared regions), 
elastic and inelastic differential cross-sections, and virial coefficients data. In parallel, these 
systems have stimulated a remarkable interest in the theoretical community, providing 
the possibility of carefully testing ab initio methods. For example, ArH2 has become a 
popular benchmark for theory, especially for calculations of rotational and vibrational 
pre-dissociation effects j2|, as for this system the widest range of spectroscopic, coUisional 
and bulk property data are available j^. 

In recent years the advent of Bose-Einstein condensates of atoms and molecules has 



revitalisec. 



physics ares rapidly evolving and attracting a considerable interest in the scientific 
community. This is due to the possibility of cooling atoms and molecules at temperatures 
below 1 K, which creates the capability of experimentally controlling chemical reactions, 
allowing accurate studies of ro-vibrational energy transfers and coherent chemistry. A 
rather promising experimental route being pursued is the creation of ultra-cold molecules 
by means of sympathetic cooling with other ultracold species. Consequently, a number of 
groups have begun the investigation of ultra-cold collisions and chemical reactions from 
the theoretical point of view [5;!], calculating ab initio elastic and inelastic cross-section, 
branching ratios, resonance energies and widths, etc. In particular, the calculation of cross 
sections (both elastic and inelastic) between different species is hugely importance for the 
planning and realisation of experiments with ultra-cold species. 

The main ingredient in an ab initio calculation is the Potential Energy Surface (PES) 
modelling the interaction among the atoms, within the Born-Oppenheimer approximation. 
In recent years, a large number of PESs have been produced for the five stable Rg-H2 
complexes (Kg = He, Ne, Ar, Kr, and Xe). Initially, the most accurate PESs were produced 



the interest in van der Waals complexes. The areas of cold and ultra-cold 
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semi-empirically, by means of fits to measured data jl| . With the advances in computational 
power and quantum chemistry methods, purely ab initio PESs of comparable accuracy are 
now available. In particular, extensive work has been done on Ar-H2 , , which is 

one of the most thoroughly studied atom-diatom complexes. More recently, the He-H2 PES 
has attracted a considerable attention. Due to the weakness of the He-H2 interaction, the 
existence of a bound state for this complex has long been debated. Recently, Kalinin et 
al. have produced evidence for its existence by means of a transmission grating diffraction 
expsHment gin parallel, a fully a, .nU,o study of its He-H, PES has been carded with 
high accuracy ^0|, which has inspired a number of theoretical papers on the vibrational 
spectrum of the complex. 



Despite the very large availability of PESs, relatively few groups have reported accurate 
quantum mechanical calculations for Rg -H^ complexes. The He-H2 system has been object 



of three theoretical papers 
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131] aiming to assess the properties of its halo state. 



following Kalinin et a/.'s experiment , whereas He-H2 collision processes at low temperature 
were studied by other groups jl^ . [isl . The ro- vibrational spectrum of Ar2 was calculated 



variationally by Moszynski and co-workers 
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were investigated by Uudus and co-workers 



, whereas H2 and D2 collisions with Ar atoms 



ITj . The rotational- vibrational spectrum of 



Kr-H2 was analysed by Wei and co-workers P[and by Zhou and Xie [l9| . The former have 
produced a highly accurate semi-empirical PES, also performing all dimensional calculations 
of the molecule's vibrational-rotational levels, whereas the latter reported a purely ab 
initio PES calculated at the CCSD(T) level of theory together with the calculation of the 
associated vibrational-rotational levels. 



The object of this work is a comparative analysis of Rg-H2 complexes, using, where 
possible, more than a single PES for each molecule. Given the very large number of PESs 
available in literature, our aim is not to provide a fully comprehensive study of them all. 



We chose to selecting some of the most popular semi-empirica' 



as those of Ro dwell and S coles 



20], and Le Roy and Hutson 



PESs from the 1980s, such 



211 ]. and all the most recent 



ones, including the one by Boothroyd et al. lOj for He-H2, Bissonnette et al. [3] for Ar-H2, 
and a few others. Although some of the surfaces are known in literature with proper names, 
each of them were assigned a new acronym in this paper in order to uniform the notation 
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among the different complexes. The molecules are studied by variational methods, by using 
the Rayleigh-Ritz and Kohn principle for bound and continuum states, respectively. The 
trial wavefunction used is based on a simple, but effective model, which assumes that the 
Rg atom rotates/vibrates around a the H2 molecule stiff in its vibrational-rotational ground 
state. This model has been used in variational many-body calculations [22], and Gianturco 
et al. proved it sound also in comparison to very accurate full dimensional calculations 
for the halo He-H2 state [l2^. We extend their work to also study cold scattering, with 
applications to the five different Rg-H2 complexes. 



This paper is organised as follows. The next Section will briefly present the model used 
to study the Rg-H2 systems, with details of the computational procedure. In Section III 
the results obtained for the five complexes are presented and discussed individually. The 
last Section is devoted to the comparative analysis of the different complexes, as well as the 
conclusions. This Section presents the elastic cross-sections at cold energy, as well as the 
structural and energetic properties of the vibrational molecular states. 



II. COMPUTATIONAL DETAILS 



The H2 molecule is strongly bound compared to RgH2 binding. Using the H-H potential 



from Ref. 23|], the lowest rotational levels in the ground vibrational band result to be 118.4 
cm~^, 354.2 cm~^ and 705.1 cm~^ for the states with £ = 1, £ = 2 and £ = 3, respectively, 
where the energy of the ground state (£ = 0) is taken as zero energy point. In this type of 
problems, it is common practice to write the system wavefunction in terms of the diatom 
molecular states 24] 

(1) 



LM ' 



^ = E 0.,(r)/.,(i?) [Y,{f)Y,{R) 

where (puei'f') is the radial part of the H2 wavefunction with u and i its vibrational and rota- 
tional quantum numbers, respectively, Y({r) and Yrj{R) are spherical harmonics associated 
to the Jacobi vectors r and R, respectively, and fur]{R) is a set of unknown functions to be 
determined. The molecule is considered in a state of total angular momentum L, and z-axis 
component M. Due to the weakness of the Rg-H2 potential we can restrict our analysis to 
L = 0. The heavier Rg atoms form compounds with L > but we will ignore them for 
the purposes of this work. For the particular case of L = we have that £ = rj and the 



dependence of \1/ on the four polar angles f and R can be reduced to just one non-trivial 
coordinate by integrating out the Euler angles 

Y,{r)Y,{R)\^ = P,{^^) (2) 

where = (r) ■ (i?) and is a Legendre polynomial. 

At this point it is possible to notice that the expansion of eq. [T] involves H2 states with 
high internal energy if compared to the Rg-H2 interaction. In fact, the vibrational energies 
of the Rg-H2 complexes are 10-100 times smaller (vide infra). It is reasonable to expect that 
the presence of the Rg atom will only slightly affect the H2 diatom, even a short distances. 
The Rg-H2 complex can thus be reasonably modelled, at low temperatures, by assuming that 
the Rg-Il2 potential only slightly perturbs the H-H Hamiltonian. Under this assumption, 
the two II2 quantum numbers v and I can be considered good quantum numbers for the 
three-body system. In mathematical terms, this corresponds to retaining only the first term 
of the expansion of eq. [H A presumably very accurate trial wavefunction for those complexes 
can thus be written as 

vI/ = 0oo(r)/(i?)/i?Po(/^), (3) 
where f{R) is an unknown function to be determined. 

A very close approach was already shown to be surprisingly accurate by Gianturco et al. 



12], with an application to the binding energy of the HeH2 system. In this work, we are 
interested in studying how well the five RgH2 complexes can be modelled by using the trial 
wavefunction of eq. ([3]). We are interested in calculating not only the vibrational energies, 
but also the the associated wavefunctions and their structural properties. Furthermore, we 
can apply the same trial wavefunction to study the continuum energy region, in order to 
calculate the atom-dimer scattering length, the effective range, the phase-shift over a range 
of energies, and finally the Rg-H2 cross sections. Those calculations will finally allow us to 
make a comparative analysis of the scattering properties in the five different complexes and 
of the several available PESs. 



A general expression for the Rg-H2 intermolecular PES is the following: 

V{r,R,^l) = V^^_^,{r,R,^l) + Vmi{r), 
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(4) 



where Vhh is the H-H potential when the Rg atom sits at an infinite distance from H2, 
and l^Rg-H2 parameterises the weak interaction between H2 and the Rg atom. The oldest 
potentials used in this work assumed a stiff H2, and are expressed in terms of R and n only jl| . 



In eq. (12]) the only unknown is the function f{R) associated to the motion of the Rg 
atom with respect to H2. If the variational principle is applied to the trial wavefunction 
one is led to solving the one-dimensional equation 



2/x dR^ 



+ Veff{R)-E 



f{R) = 0, (5) 



where 



Veff = J (PUr)P^{li)V{r, R, fi)drdfi, (6) 



and E is the energy of the three-body system with respect to the H2 binding energy. In par- 
ticular, the zero point for the energy axis is chosen with the Rg atom at an infinite distance 
from the II2 molecule in its vibrational-rotational ground state. All energies reported will use 
this convention. We will also use cm~^ as units for energies and Afor distances. Equation ([5]) 
is solved using the Numerov algorithm, and the effective potential is calculated on a uniform 
grid of 10,000 points equally distributed between 0.01 and 100 A. The equation above can 
be solved for negative values of E, yielding upper bounds to the vibrational energies, and 
for positive values of E, yielding the complex scattering wavefunction Je- In this latter case, 
a first order estimate for the observable of interest, say C, is obtained from the asymptotic 
part of the function fE{R), and the second order estimate by means of the Kohn variational 
principle: 

^2nd ^ ^ist _ ^^^1 ^ y^^^^j^^ _ ^ |^^^_ (7) 

With the procedure just described the phase-shift can be obtained over a range of 
energies E. The centre of mass elastic cross section is then expressed as 

, ^, , sin^ b , , 

a{E) = 47r-^, (8) 



where k = \j2ij,E/h . At low energy the effective range expansion for 6 reads [25[ 

A; " a, 2 ' 
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where is the scattering length and Veff the effective range. Consequently, the limit for 
elastic cross section at zero energy is 



a{E = 0) = Awal- 



(10) 



The effective range re// is obtained by means of the integral relation 25 1 

2 



fo{Rr-{R-as 



(11) 



where /o(-R) is the solution of eq. ^ ioi E = 0, with the appropriate normalisation. 



In order to test whether the ansatz in eq. ([3]) for the wavefunction is a good choice, 
it is possible to calculate the contribution of the next H2 channel in energy to the desired 
observables. Namely, a two channel trial wavefunction which includes the H2 {u = 0,i = 2) 
state can be build as follows: 

vl> = 0oo(r)/(i?)/i?Po(/^) + (Po2{r)g{R)/RP2{fi), (12) 

where / and g are two unknown functions to be determined . This choice is motivated by 
energy and symmetry considerations, as the state considered is the allowed H2 rotational- 
vibrational state closest in energy to the the ground state, as the state with i = 1 belongs 
to a different symmetry block. The implementation of the variational principle leads to a 
system of two one-dimensional coupled differential equations. In order to determine / and g, 
they are expanded in a polynomial basis, and the coefficients of the expansion are determined 
by solving a generalised eigenvalue problem for the bound state or a linear problem for the 
scattering states. Results associated with the trial wavefunction [12] will only be presented 
for a specific case, in order to show that the contribution of the state (z/ = 0,£ = 2) is very 
small. 



III. RESULTS 



A. He-Ho 



The He-H2 system is the computationally most interesting one due to the presence of a 



halo state. Halo states 



are very weakly bound and extremely sensitive to the minimum 
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details of the PES. As the He-H2 van der Waals interaction is a consequence of big cancella- 
tions between the different terms of the molecular Hamiltonian, very accurate calculations 
of the PES are needed in order to provide a reliable ab initio PES and a subsequent estimate 
of the HeH2 binding energy. The PES of Bissonnette et al [10] (heh2c) was computed fitting 
a very large number of ah initio energies calculated at the MR-CISD (multiple-reference 
single and double excitation configuration interaction) level of theory. Other available PESs 
for this complex are the semi-empirical from Ref. j20| (heh2s), and from Ref. 27] (heh2b). 
Table [T] shows calculated observable for the HeH2 complex, namely the parameters Rmin, 
e and Rq, representing the position of the minimum and the well depth of the effective 
potential K//, and the point where it crosses the i?— axis, respectively. The HeH2 system 
has only one bound state, of energy E. Other listed quantities in Table [I] are the mean 
values of the kinetic and potential energy, (K) and {V), respectively, and the mean values 
of the Jacobi coordinate R, for the ground state wavefunction. Scattering parameters listed 
are the scattering length a^, the effective range rg//, and the value of the cross-section at 
zero energy a{E = 0). As in other works, we have made calculations, only for PES heh2c, 
for both isotopomers '^HeH2 and ^HeH2. The last two rows of Table [T] show the correction 
to the ground state energy and to the scattering length obtained by using the two-channel 
wavefunction of eq. (fT2!) . 

The three PESs used show good agreement with each other. The biggest difference, as 
expected, is in the energy of the bound state, where differences are up to a maximum of 
20 %, whereas differences in the scattering lenght are less than 5%. The only vibrational 
bound state supported shows a highly correlated structure, reflecting its halo character. Its 
binding energy is given by a big cancellation between the kinetic and potential energies, 
and both its size and the scattering length are larger than the range of the He-H2 potential. 
Those characteristics are conflrmed by all three surfaces considered. A similar structure, yet 
further enhanced, is present in ^HeH2, which, being lighter, is less bound. The last column 
of Table [T] displays results for ^HeH2 calculated with the heh2c PES. 

Table [IT] compares the results obtained in this work with the literature. For both •^HeH2 
and ^HeHo there is a very good agreement with the more accurate calculations of Xiao and 
Poireir ISj, and of Gianturco et al 12|, relative to the heh2c PES. From the latter it is 
shown both the value obtained with a three-dimensional Discrete Variable Representation 
expansion and the value (indicated as He-2H ) obtained with an approach similar to the one 
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used in this work. Results from Balakrishnan et al. IJ], relative to the PES heh2b, also 
compare quite well with the ones from this work for both He isotopes. 
For this complex we have also calculated its binding energy and scattering length, relative 
to PES heh2c, using the more elaborate trial wavefunction of eq. f[T^ . We obtain the values 
of —0.0369 cm~^ for E and 22.648 for a^, which allow us to conclude that the contribution 
of neglected terms in the ansatz function ([3]) is negligible. 



B. Ar-H, 



The attention devoted to this complex in the past years by the quantum chemistry com- 
munity has resulted in the availability of a very large number of PESs. The most commonly 
used is the semi-empirical PES of Ref. Q (arh2f). We have also used the surface of Scoles 



20[ | (arh2s), Le Roy and Hutson (Ref. 2l|], arh2b), and Fit 2 of Schwenke et 



and Ro dwell 

al (arh2c). Table UTTl lists all relevant parameters for this system. Due to the depth of 
the effective potential, two vibrational states are present (n = 0, 1). The ground vibrational 
state (n = 0) displays a strongly bound character, as for example its binding energy is of 
the same order of magnitude than its kinetic and potential energies, and the variance of R is 
much smaller than its mean value. Conversely, the excited vibrational state n = 1 displays 
a slightly weaker character, which will be analysed graphically in Fig. [H altough it can not 
be considered a halo state. Uudus et al [17] report quantum mechanical calculations of cross 
sections for ArH2 using the arh2c and arh2f PESs. In particular, their calculated scattering 
length, which is 9.95 Afor arh2c and 10.09 Afor arh2f agree well with the ones in Table UTTl 
All PESs agree within a factor of 10%. The inclusion of the second channel in the trial 
wavefunction brings a small contribution to both the binding energies and the scattering 
length, of less than 1 %. 



C. Ne-H2, Kr-H2 and Xe-Hs 



Ref. 



(neh2b), from 



For Ne-H2 there are no recent PESs. We used the surfaces from Ref. p 

29l | (neh2c), and from Ref. 20|| (neh2s). Results for this complex are shown in the 



first three columns of Table IIVI The three surfaces agree within 20%. They support a 
single vibrational band, relatively deep in energy (about 20%) of the potential well. The 
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halo structure present in HeH2 is lost in this system. Tennyson and Sutcliffe |30!] using PES 
neh2b obtained an energy of 4.67 cm~^ for the vibrational band origin, in agreement with the 
result of this work. The calculation of accurate surfaces for the heavier atoms, Kr and Xe, 
is made more difficult by the rapid increase of the number of electrons. For the latter, the 
only available surface was determined by Le Roy and Hutson 21j (xeh2b). The associated 
results are shown in the last column of Table HVl For Kr, we used a semi-empirical PES from 
;he same authors (krh2b), a much more recent semi-empirical PES by Wei and co-workers 
isl (krh2c), and the fully ab initio surface of Zhou and Xie 19] (krh2d). Wei et al report 
calculated values for the two band origins, of 28.468 cm~^ and 1.653 cm^^, which agree with 
the ones presented in Table IIVI Also Zhou and Xie [3] report their calculated values for 
the two band origins, which are 25.688 and 1.048, again in good agreement with the ones 
calculated in this work and reported in Table IIVI The agreement between PES neh2b and 
neh2c is rather surprisingly good, with difference of less than 5% also in quantities, like the 
scattering length, which are sensitive to every detail of the surface. This agreement leads us 
to believe that also PES xeh2b would prove reliable if compared to more accurate surfaces. 
On the other hand, differences between those two PESs (krh2b, krh2c) and krh2d are quite 
large, and are a consequence of the more shallow depth of the latter. In both Krll2 and 
XeH2 the two bound states exhibit a deeply bound structure, shown by the mean values 
of the observables presented in the Table. In fact, the binding energy is comparable to the 
kinetic and potential energies, the associated radial distrbution is rather peacked around its 
maximum, and the scattering length is rather small. For all complexes the inclusion of the 
second channel changes the results of an insignificant amount if compared to the uncertainty 
associated with the PES. 



IV. COMPARATIVE ANALYSIS AND CONCLUSIONS 
A. Comparative analysis 

Figure [T] shows the radial distribution p{R) of the Rg atom with respect to the H2 centre 
of mass, for selected RgH2 complexes and vibrational states. In particular, the halo states 
associated with the two IIeH2 isotopomers, and the two vibrational states of Arll2 have 
been plotted. Note the bi-log scale of the plot. The distribution densities associated with 
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the remaining RgH2 complexes display characteristics very similar to the ground vibrational 
state of ArH2, and are therefore not shown. The ground state of ArH2 displays a rather 
peaked structure, and the function p{R) is different from zero only in a small range of values 
of R near its peak. Conversely, the density of the two HeH2 isotopomers show a completely 
different structure, being more flat and extending over a broad range of values of R, well 
beyond the range of the effective He-H2 potential. Those are typical features of a halo 
state. The halo structure is more pronounced in ^HeH2 than in ^HeH2, as expected by the 
ligMe. .a. of the fo™. Q. The excited vibrational =tate of AtH. shows a sttuctu. 
midway between the two just illustrated: it is quite broad in R but less than the HeH2 
complexes. At short distances, all distributions fall rapidly to zero as a consequence of the 
hard repulsive core of the inter-molecular potential, which is consequence of the closed-shell 
character of its constituents. At zero incident energy, the largest elastic cross sections are 
found in He-H2 and Ar-H2, with the others becoming smaller the heavier the Rg atom. 
This behaviour can be easily interpreted following the analysis of the radial distributions 
of the bound states, as the two complexes have a halo bound state. It is well known that 
for halo complexes the scattering length is inversely proportional to the square root of the 
binding energy (a^ ~ 1/V^), therefore the smaller the latter the bigger the former. Elastic 
cross sections for a range of energies are shown in Fig. [2l In the limit of zero energy, the 
difference between the smallest (Xe-H2) and the largest (He-H2) is more than two orders 
of magnitude. However, in commenting those results one should always recall that a small 
change in the PES can result in a very large modification of the cross section for the energy 
range considered, due to the extreme sensitiveness of the scattering length to the strength 
of the interaction. In order to have a rough check of the sensitivity of each complex to small 
modifications of the PES, we have arbitrarily modified each PES by a factor 1 + A, 

VxiR) = il + X)V,ffiR), (13) 

where A = recovers the real case. This same procedure was applied to test the Efimov 



character of the excited state of the Helium trimer 



31 



32| . Figure [3] shows the changes in the 
Rg-H2 scattering length as a function of A, compared with the absolute value of the ^He-H2 
scattering length. The Ne-H2 and Kr-H2 complexes display very little sensitivity to changes 
in the PES of up to 25%, in both directions (weakiniing and strengthening), whereas Ar-H2 
and Xe-H2 are much more affected by changes to their PESs. In particular, the former will 
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present the same zero energy cross section of ^He-H2 when its PES is weakened by a factor 
1/5, whereas the latter when its PES is strengthened by a similar amount. In percentage 
terms, this implies that a modification of up to 20% of the Xe-H2 PES can change the 
associated zero energy cross section of up to two orders of magnitude. 



B. Conclusions 



In this paper we have studied the vibrational band origins and ultra-low energy scattering 
properties of the Rg-H2 complexes, Rg=He, Ne, Ar, Kr and Xe. A simple one-channel trial 
wavefunction based on a independent H2 model has been used to investigate both the discrete 
and cold energy range, in conjunction with variational methods. A number of PESs were 
used to mimic the Rg-H2 inter-molecular interaction. The accuracy of this approach has been 
checked in a number of ways. Firstly, the contribution of the next H2 channel was evaluated 
by implementing a two-channel calculation. This larger calculation was performed for all 
PESs to calculate the binding energies and scattering lengths for all complexes considered. 
It showed that the inclusion of the second channel brings a very little contribution, of less 
than 1 %, which is insignificant if compared to the uncertainty associated to the molecular 
PES. Secondly, the results obtained compare well with other published works where more 
accurate full-dimensional approaches were pursued. The different PESs agree reasonably 
with each other, in particular the semi-empirical sets of Rodwell and Scoles [2^ and Le 
Roy and Hutson 2l(] yield results close to the more recent, and presumably more accurate, 
potentials. The Rg-H2 elastic cross section at ultra-low energy varies by up to two orders of 
magnitude when the different Rg atoms are considered. He-H2 has the greatest cross section 
due to the presence of the halo state in the trimer, whereas Xe-H2 has the smallest one. The 
Ar-H2 elastic cross section is also large, possibly due to the presence of the relatively weakly 
bound excited vibrational band in this complex. Finally, we have checked the dependence of 
the atom-dimer scattering length on the strength of the atom-diatom effective potential. It 
can be argued that the uncertainty in the Rg-H2 PESs can qualitatively affect the conclusions 
of this study, as the zero-energy scattering length can be dramatically influenced by small 
(or even by tiny) adjustments of the intermolecular PES. To test those effects, the Rg-H2 
PES has been arbitrarily modified by a factor (1 + A), and the effect of varying A on the 
atom-dimer scattering length has been studied. Whereas Ne-H2 and Kr-H2 have shown 
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stability over the range considered (—0.25 > A < 0.25), Ar-H2 and Xe-H2 proved very 
sensitive to such changes, and, for instance, the Xe-H2 cross-section increases by two orders 
of magnitude if the PES is strengthened by 20%. As the Xe-H2 PES is also the least well 
known among the five complexes considered, this result suggests that a more detailed study 
of the Xe-H2 interaction is needed. 

In conclusion, the proposed model can be easily extended to the study of four-body ultra- 
cold and cold collisions of the type Rg-Rg-H2. In fact, also in this case the H2 molecule can 
be approximate as a single particle, and the importance of inelastic collisions leading to the 
formation of Rg-H2 complexes in the four-body collisions can be studied by means of three- 
body techniques. Regarding applications to Rg-H2 in higher angular momentum states, or to 
different complexes of experimental interest, such as Rg-NHs, the simple trial wavefunction 
proposed in this work may be not as effective as in the particular case considered here. Even 
for those systems, however, it will be interesting to compare results calculated with a many 
channel and a single channel calculation, recalling that the largest source of error is usually 
given by the uncertainty in the PES. 
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heh2s[20] heh2b[27] 


heh2c[10l 


Rq 


9 QQA 




2.993 


-RfTliTl 


(J. (J 1 tJ 


3 4?0 


3.373 


e 




8 450 


9.499 




4TTprT 








E 


-U.Uooo 


-U.UouU 


-u.uouo 




(K) 




49fiQ 


0.4831 


0.1564 








-0.5198 


-0.1599 


{R) 


ZD. -LO'4: 


1 A Am 


13.372 


32.922 


V \ / 


32.076 


17.711 


16.306 


41.586 


Us 


23.647 


24.708 


22.690 


67.638 


^eff 


5.08 


5.16 


5.00 


6.08 


a{E = 0) [A2] 


7027 


7672 


6470 


57500 


AE 


-0.0002 


-0.0002 


-0.0001 


-0.00005 


Aas 


-0.069 


-0.085 


-0.052 


-0.478 



TABLE I: List of properties of the HeH2 complex. All distances are in A and energies in cm~^. 
For all potentials the results presented refer to the isotopomer ^HeH2, whereas for the potential 
heh2c the results are reported for ^HeH2 in the left column and for '^HeH2 in the right column. 
The upper part of the table reports quantities associated to the He-H2 effective potential, namely 
the crossing point with the i?— axis Rq^ the minimum Rmin and the well depth e. The upper 
middle part reports the energy of the bound state E, and the mean values of the kinetic energy K, 
the potential energy V , and the intermolecular distance R, along with its square R?. The lower 
middle part reports the scattering length a^, the effective range rg// and the zero energy cross 
section a{E = 0). Finally, the bottom part reports the correction to the binding energy and to the 
scattering length obtained by using the two-channel trial wavefunction with respect to the single 
channel one {Ax = x^^-* — x^^-*). The values used of the atom's masses are 1.00794 amu, 4.00260 
amu and 3.01603 amu for H, ^He and ^He, respectively. 
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heh2c[10] 




heh2b[27] 




Ref. [13j 


Ref. [12] 


this work 


Ref. [14] this work 






He-Hs He-2H 






^He-H2 


-0.03640 


-0.03634 -0.03803 


-0.0368 


-0.0298 -0.0300 


3He-H2 


-0.00327 


-0.002916 


-0.00348 


-0.0016 -0.0019 



TABLE II: Compilation of the values for the binding energies, in cm^^, for ^HeH2 and ^HeH2 
present in the literature, compared to the values obtained in this work. The agreement is good for 
both isotopomers and for the two PESs considered. 
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arh2s[20] arh2bf21] arh2c[7j arh2f[341 


Rq 


3.144 


3.179 


3.194 


3.179 


Rmin 


3.551 


3.580 


3.610 


3.591 
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53.922 


50.88 


50.982 


50.479 




n=0 
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-23.952 


-22.113 - 


-22.205 


-21.975 
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11.156 
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4.078 
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3.142 
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(^) 
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-3.621 


-3.205 
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7.752 
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7.940 


\/(^') 

V \ / 


7.680 


8.140 


8.466 


8.364 


as 


8.710 


9.523 


10.102 


9.930 




6.12 


5.22 


5.22 


5.21 


a{E = 0) [A2] 


953 


1140 


1282 


1239 


Ef^ 


-0.016 


-0.016 


-0.015 


-0.015 




-0.007 


-0.002 


-0.001 


-0.001 




-0.002 


-0.009 


-0.010 


-0.009 



TABLE III: List of properties of the ArH2 complex. All distances are in A and energies in cm~^. 
The mass of Ar has been set as 39.948 amu. Results are shown for the ground (n=0) and excited 
(n=l) vibrational states. For the explanation of the different reported quantities see caption of 
Table H 
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Ne-H2 






Kr-H2 




Xe-H2 




neh2b[28l 


neh2c[29j 


neh2sl"20] 
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xeh2b[21j 
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2.963 
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3.304 


3.331 


3.510 




3.300 


3.342 


3.300 


3.720 


3.730 


3.759 


3.940 


e 


22.978 


21.071 


23.565 


58.849 


58.466 


54.517 


64.914 




n=0 


E 


-4.741 


-3.959 


-4.851 


-28.453 


-28.493 


-25.648 


-33.194 


(K) 


5.448 


4.917 


5.563 


12.203 


12.090 


11.518 


12.949 


iv) 


-10.189 


-8.876 


10.414 


-40.657 


-40.583 


-37.166 


-46.142 
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V \ / 
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n=l 


E 








-1.734 
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-3.139 
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8.996 
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(R) 
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6.434 


6.890 


6.153 


V \ / 








6.585 


6.614 


7.120 


6.278 


as 


3.360 


3.846 


3.309 


5.509 


5.659 


6.961 


1.819 


Tpff 
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14.37 


10.70 


14.92 


14.83 


13.94 


9.65 


288.78 


a{E = 0) [A2] 
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138 


381 


402 


609 


41.6 




-0.006 


-0.006 


-0.007 


-0.024 


-0.021 


-0.024 


-0.030 










-0.003 


-0.002 


-0.003 


-0.005 




-0.004 


-0.003 


-0.003 


-0.006 


-0.006 


-0.007 


-0.012 



TABLE IV: List of properties of the XH2 complexes, with X= Ne, Kr, Xe. All distances are in 
A and energies in cm^^. The values used for the atom's masses are 20.1797 amu, 83.8 amu and 
131.293 amu for Ne, Kr and Xe, respectively. Results are shown for the ground (n=0) and excited 
(n=l) vibrational state. For the explanation of the different reported quantities see caption of 
Table H 
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Figure Captions 
Figure 1. 

Density distribution P{R) for selected RgH2 vibrational states. In particular, the density 
distributions for the bound states of ^HeH2 and "^HeH2, and the two vibrational states for 
ArH2 are plotted. The densities shown refer to calculations made using the heh2c and arh2b 
potentials. Note the bi-log scale. 
Figure 2. 

Elastic cross sections for the five Rg-H2 complexes, plus the ^He-H2 isotopomer. The 
energy is given with respect to the centre of mass frame, ranging from 10^^ to 10 cm^^. 
Note the bi-log scale. The curves shown are: dotted line, dot-dash line, dot-dot-dash line, 
dash-dash-dot line, dashed line and continuum line for ^He-H2, 43He-H2, Ar-H2, Kr-H2, Ne- 
H2 and Xe-H2, respectively. The displayed cross sections are relative to phase-shifts obtained 
with the heh2c, neh2s, arh2b, krh2b and xeh2b potentials. 
Figure 3. 

Dependence of the scattering length a^ on the parameter A (eq. 0131) ). The two horizontal 
continuum lines indicate the threshold where a^ reaches the same value in absolute terms 
of the ^He-H2 scattering length. Calculations have been made using the same potentials as 
indicated in the caption of Fig. [2l 
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